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1^ ■ Abstract 

, A particle-hole model is developed to describe the excitation spectrum of single A hypernuclei and the 

possible presence of collective effects is explored by making a comparison with the mean- field calculations. 
Results for the spectra of , ]yO , j^Ca , ^°Zr and j^^Ph hypernuclei are shown. The comparison 
with the available experimental data is satisfactory. We find that collective phenomena are much less 
important in hypernuclei than in ordinary nuclei. 

PACS: 21.80.-t-a 



1 Introduction 

OO ' The properties of hypernuclei have been widely studied in several experiments and today they are the 

, object of a rich experimental program. The first experiments were aimed at identifying the hypernuclei 

0^ ' and at determining the energy of the hyperon embedded in the nucleus. Presently, the improvement of 

, the experimental techniques allows the measurement of high-quality excitation spectra, which have been 

^ ■ studied in {11+ ,K+) [Dllllll], {K-,tt-) [51E1[7] or {e,e'K+) ^M\M reactions. In the future, further 

advance in these studies will be gained by means of 7-ray spectroscopy, which will become a standard 
investigation tool [Til US] • 

From the theoretical point of view, single-A hypernuclei have been studied within the mean-field (MF) 
;h ' approximation: the hypernucleus is seen as a many-body system containing an impurity, and its properties 

5^ , are solely determined by those of the hyperon moving in an average potential generated by the interaction 

with the nucleons. When applied to various medium heavy nuclei, the MF model is quite successful in 
reproducing the single hyperon energies, which are single particle properties of the system. On the other 
hand, the excitation spectra of hypernuclei are, in general, related to the system as a whole. 

In purely nucleonic nuclei there are large differences between the experimental excitation spectrum 
and that predicted by the MF model, where the excitation energies are simply given by the differences 
between the energies of the particle and hole states. This indicates that the many body-terms of the 
hamiltonian neglected in the MF approximation play a relevant role. It is therefore interesting to study 
also the excitation spectra of hypernuclei by using models beyond the simple MF approach. 

In the past this kind of studies have been done within the Tamm DancofF Approximation (TDA) 
I14j , which describes the many-body excited state as a linear combination of particle- hole excitations 
of a MF ground state. The theoretical inconsistencies of the TDA theory are overcome by the Random 
Phase Approximation (RPA) theory which considers some ground state correlations jTS] . 
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We have developed an RPA-like model which describes the spectra of single Ahypernuclei. In our 
model, the hypernucleus is formed when a nucleon, from its state below the Fermi surface, a hole state, 
is transformed into a A which lies in one specific particle state. The excited state of the hypernucleus is 
a combination of various particle- hole excitations. This combination is ruled by the residual A -nucleon 
interaction. The relevance of the residual interaction, and therefore the possible presence of collective 
phenomena, is investigated by comparing MF results with those obtained by our model. 

In the next section we present the details of the model. In Sect. [3] we describe how we have chosen the 
input of our calculations, namely the set of single particle wave functions and energies, and the effective 
A -nucleon interactions. A selected set of results concerning hypernuclei obtained from doubly-closed shell 
nuclei is discussed in Sect. [H Finally, in Sect. [5] we summarize the basic findings of our study and draw 
our conclusions. 



2 The formalism 



The starting point of our model is the MF description of the many-body system. Each baryon is described 
by its single particle wave function 0^ and energy e^, where a indicates the set of quantum numbers 
which identify the single particle level. 

The definition of the single particle basis allows us to write the many-body hamiltonian as [15] 



-ff = e«ala„ - ^ X! ^hh'hh' + | X! ^aa'P/S'N 



al,a'l,ai3>ap 



(1) 



hh' 



where we have indicated with and a the creation and annihilation operators, and with N the normal 
ordered product operator |16j . In the above equation, the bar over the interaction V indicates that both 
direct and exchange terms are considered, the h label identifies a state below the Fermi surface, while 
the Greek subindexes refer to states above or below the Fermi surface. 

We describe the formation of the hypernucleus as a time dependent fluctuation of the MF ground 
state. We solve the time dependent Schrodinger equation by using the variational principle 



S\<^{t)\ 



H -ih 



|*(t)> =0 



(2) 



We search for the minimum of the above energy functional within the Hilbert subspace spanned by 
trial wave functions of the form 



_ -lEot/h 



exp 



Ah 



1*0 > , 



(3) 



where |<I>o > indicates the MF ground state of the system, formed by nucleons only. In equation ([3]) we 
have indicated with aj^ the operator which creates a A , and with ah the operator annihilating a nucleon. 
The coefficients Cf^h{t) are complex numbers. 

By using the ansatz ([3]) in Eq. ([2]) we express the variational equation as 



d 



dClhit) 



^-4 



>= 



(4) 



By making a power expansion of the exponential of Eq. ([3]), and retaining only the terms up to the 
second order in C{t), we obtain the expression [T5] 



CAh{t){eA — ih) + CX' h' it)V AA' hh' 

A'h' 



J2CA'h'it)VAh' 



hA' 



A'h' 



(5) 
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Since we are looking for oscillating solutions, we write the coefficients as 

C^f,{t)^W^he'"'' + ZX„e'''' , (6) 
and by equating positive and negative frequencies components we obtain tlie equations 

WAh{eA - e/i) + ^ VAh'hA'WA'h' + ^ VAA'hh'ZA'h' = hiyWA'h' (7) 

A'h' A'h' 

ZhX^K - + E y^h'hA'Zl,j,, + Vaa'hh'WI,,,, - -hvZl,,, . (8) 

A'h' A'h' 

Introducing the matrix elements 

HAhA'h' = i^A — £h)5AA'Shh' +VAh'hA' (9) 
lAhA'h' — VAA'hh' , (10) 

equations ([7]) and ([8]) can be expressed in the following matrix form 

where we have defined u = hv. We have rewritten the secular equations in a matrix form to show their 
similarity with the traditional RPA secular equations J7j . Our derivation of the secular equations shows 
that the W and Z amplitudes are introduced on equal footing, and that the choice of neglecting one of 
the two terms of Eq. Q is a further approximation of the theory. On the other hand, the derivation of 
the secular equations done with the equation of motion method [151 117j indicates that the Z amplitudes 
are related to ground state correlations. Neglecting the Z amplitudes produces TDA secular equations 
which have been used in the past [131 HI] to describe hypernuclei. A discussion on the relevance of the 
ground state correlations introduced by the Z amplitudes will be done at the beginning of section [H 
However, we would like to point out here that the ground state correlations we include are only those 
described by A-nucleon excitations, and not by nucleon-nucleon excitations. For this reason our model 
does not include nucleonic correlations as it is done in the traditional RPA theory. 

In our calculations the ground state of the nucleus and the final state of the hypernucleus have definite 
angular momentum values. For this reason we rewrite the above equations in angular momentum coupling 
scheme. The angular momentum coupled amplitudes are defined in terms of the uncoupled ones as 

Wih = < UmAjhmhlJM >WAh (12) 

'niAmh 

Z'L = <UmAjhmh\J-M>{~lY-^'ZAh. (13) 

where the j's indicate the total angular momentum characterizing the single particle wave functions and 
the m's their z-axis components. 

We express the terms ([9] - ITO)) of the secular matrix (|TT|) as 

H'lhA'h' = {<^A - Qi)6aA' Shh' +VAhA'h' (14) 
I Ah A'h' — Ah A'h' J (15) 

where we have defined 

,,J ^.Jdir ^,Jexch 



J _ Jdir Jexch 

^AhA'h' — ^ Ah A'h' ~ ^AhA'h' 

\^(_l)Jh+iA'+-ff^ i jh -J 

I J A' jh' I'' 

[< jAjh'K\\V\\jHjA'K > -(-ly^+lA'-K < j^j^,K\\V\\jA'jhK >] , (16) 
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and 

UAhA'h' = ^AhA'h' - UAhA'h- = (-1) ^^A/ih'A' ■ (I'j 

Here we use the notation 

J = \/2J+ 1 , (18) 

and the term in curly brackets is the Wigner 6j coefficient. The double bar in the matrix elements in 
Eq. ()16p indicates that in the calculation of the angular part we consider the reduced matrix element as 
defined by the Wigner-Eckart theorem [T5] . 

In calculating the interaction matrix elements we use the Fourier transform of the interaction expressed 
in momentum space 



(27r)3 

[FA{q) + GA{q)cT{l)-(T{2) + HA{q)Si2{ci)] , (19) 

where we have used q = |q| and we have indicated with cr the Pauli spin matrices and with 6*12 the usual 
tensor operator [18]. The above expression allows us to separate the terms dependent on ri from those 
dependent on Y2 ■ The explicit expressions of the matrix elements for the various terms of the interaction 
are given in the Appendix. The solution of Eq. ()lip provides the excitation energies lo and also the wave 
function of the excited state. 

We then consider the transition from the nucleus ground state to an excited state of the hypernucleus 
induced by a generic one-body operator Tj. The transition matrix element is given by 



< j\Tj\o >= j2 r^h- < > 

Ah 

+ {-ly+'^-^'^Zi^ < MTjWja > ] . (20) 

For the natural parity states we use the operator 

Tj=jj{qr)YjMm, (21) 

where jj{qr) are the spherical Bessel function, and Yjm the spherical harmonics. The matrix elements 
to be inserted in Eq. (PO)) are 

< JaWTjU >= J drr^Rl{r)Rk{r)3.j(qr) < ja\\YjU > , (22) 
where R{r) is the radial part of the single particle wave function and 

< ja\\Yj\\jb > = < la^ja\\Yj\\lb^jb > 



(-i)^"+^^(^l ^ \ yiL+h+j) 



(23) 



In the above equation we have used the Wigner 3-j symbol .18 , and ^(Z) = 1 if Hs even, and otherwise. 
For the unnatural parity states we use the operators 

Tl=j±i = jL{qr) [YLin) (g> rTil)]ij , (24) 

therefore the matrix elements to be inserted in Eq. ((20|l are 

<3a\\TL\\jb > = 

drr^R:ir)Rb{r)jL{qr) < jaWlYUn) ^ cT{l)f \\jb > , (25) 
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with 



<3a\\[y.J+s' 



\'U> = (-1) 



la+h+]b + 



1 jajb Xa+Xb + SJ + (is,l 



An 



ia jb J 



where we have defined 



(26) 



and s = ±1. 

We also consider the transition between two excited states of the hypernucleus induced by a generic 
one-body operator {EM). By using standard techniques [T71 [12] we obtain for the transition matrix 
elements 

< $,|(£;Af)|$', > 

= E E Wlh'^i'w < M{EM)\A' > -Zl^zi,,, < A'\(EM)\A > 



Ah A'h' 



E E -^AA' [w-a^W^a:,.' < h'\{EM)\h > -Zl^zi^, < h\{EM)\h' > 



(27) 



Ah A'h' 



We calculate the electromagnetic transition probabilities between two states of the hypernucleus by 
using the expression [20] 



4k 



if 



fic2(2J, + 1)(2L+ 1) 



\{Jf\\{EM)LU) 



(28) 



where the indexes i and / indicate, respectively, the initial and final state of the transition, k is the 
modulus of the emitted photon momentum, and (EM) l is the electromagnetic operator, of multipolarity 
L. In the above equation, as before, the double bar indicates that we have to evaluate the reduced matrix 
element of the angular dependent part. 

We find that, in the angular momentum coupling scheme, the reduced transition matrix element 
between the two excited states i.e. Eq. P7|. can be written as 



J' L J 

iA jh Ja' 



< j\\iEM)L\\J' >= ii'E E I 

Ah A'h' { 

(^iy-+^'^+-''+i^wl^wi:^,<u\\{EM)L\\u, > 

-{-l)^A'+,^'+J'zi,X'h' <JA'\\iEM)L\\3A>' 
^6,A' ( ^ ^' M 

t JA Jh Jh' J 

(-1)^^' +^-+ J wl^wC, ^, < n,\\{EM) lWj^ > 
-i-l)^^+^>^'+J'+^ZLzi:,,<j,\\iEM)r^\\j,, > 



(29) 



The expressions of the single particle matrix elements for natural parity (electric), and unnatural 
parity (magnetic) transitions, are given in Ref. [21] . 
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3 Details of the calculations 



We have applied the formahsm described in the previous section to calculate the excitation spectrum of 
single Ahypernuclei formed by the excitation of nucleonic doubly closed shell nuclei. Specifically, we have 
investigated the ^^C , , fCa , fZr and f^Ph hypernuclei. 
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Table 1: Values of the parameters of the three A-nucleon interactions used in this work. The LN5 interaction is 
of zero range type, and, in the table, we indicate the values of the JJ'^'"'^ parameters of Eq. (|32|l in MeV fm^. The 
other two interactions have been constructed as indicated by Eqs. (I33|l and (|34p . The values of the w coefficients 
for the scalar (c) and spin (cr) channels are expressed in MeV, those of the tensor channel (t) in MeV fm~^. The 
ranges of the gaussians Pi are the same for both LND and LNDE interactions and their values are expressed in 
ftn. 
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Table 2: Energies in MeV of the a^O spectrum, obtained in MF approximation, and with the three interactions 
used in our work, compared to the experimental values given in Refs. [3211331 [3T| . The energy resolution obtained 
in [321 [33] is of few keV, while that of the 2+ state is of 2 MeV [31] . The MF particle-hole components are also 
shown. 



Lr analogy to what is done in the application of the Landau-Migdal theory [55] to finite nuclear 
systems [331 Hll) the mean field is described by a phenomenological potential, in our case a spherical 
Woods-Saxon well, of the form 



U\r) 



-V* 



■f exp 
he 



exp 



{r-Rl)K 



{l + exp[(r-i?L)/«L]} 



Icr 



(30) 
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where is the pion mass and the Coulomb term V^{r), active only for protons, is that produced by a 
homogeneous charge distribution. 
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Figure 1: Momentum dependence of the A- nucleon interactions for each channel of Eq. (I19p . The LNS zero- 
range interaction is constant in momentum space. The parameters of the other two interactions are given in Tab. 

HI 



For the nucleonic part we used parameterizations chosen to reproduce the single particle energies of 
the levels close to the Fermi surface and the r.m.s. charge radii. The parameters of these potentials have 
been taken from the literature [53] for the ^^C , ^^O , "'"Ca and ^"^Pb nuclei. For the ^°Zr nucleus we 
used VoP=55.88 MeV, l^o"=48.12 MeV, Vf^ = V{1=7.7 MeV, and ag'" = afj"=0.645 fm, and i?g'"=5.695 
fm, i?f;"=5.68 fm. 

In our calculations also the A mean field potential is described by a spherical Woods-Saxon well with 
a set of parameters values very close to those given in We used V^^~ 29.0 MeV, VJ^=2.0 MeV, and 

— a^=0.6 fm for all the hypernuclei, and we changed the radii in order to reproduce the ground state 
binding energies of the various hypernuclei. For , , A^Ca , ^°Zr and ^"^Pb we used Rq = i?*g=2.50, 
2.73, 4.00, 5.60, 7.60 fm, respectively. 

In both the nucleonic and A sectors the mean-field potential is diagonalized in a truncated harmonic 
oscillator basis. This procedure discretizes automatically the continuum. Since we are interested in the 
lowest excited states a proper description of the escape width is not necessary [53]. In Sect. [3] we discuss 
results obtained for excited states dominated by particle-hole transitions where the A is in a bound state. 

Our calculations have been done with a configuration space which considers 54 single particle levels 
describing the A hyperon, up to the 5pi/2 level. This corresponds to consider 10 major harmonic oscillator 
shells. A comparison with results obtained with a smaller configuration space of only 22 states, up to the 
3pi/2 level, shows differences in the eigenvalues of about 10, 20 keV. 

As already indicated in Eq. , the residual interaction between A and the nucleons is considered 
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to be a sum of scalar, spin and tensor terms 



VAN^V^^ir) + V^^{r)a{K)-<T{N) 



(31) 



In the above equation A'' indicates the nucleon and the V^]^'* {r) terms are the Fourier transforms of the 
Fa, Ga and H\ terms of Eq. (|19p . With respect to the A -nucleon interactions commonly used in shell- 
model calculations [H], we neglect the spin-orbit term. The contribution of this term is rather small [27] . 
and we consider it at the mean-field level by including a spin-orbit term in the Woods-Saxon potentials. 
The effects of the spin-orbit terms of the residual interaction have been studied in self-consistent RPA 
calculations for ordinary nuclei [28j and they have been found to be rather small. 

The sensitivity of our results to the A -nucleon interaction has been tested by comparing the results 
obtained with three different parameterizations of the V^'^'^{r) functions. In a first parameterization, 
which we call LNJ, the three functions are assumed to be of zero-range type 

VKN^r) = UrM'5{r) . (32) 

The values of the constants C/'^''^'* are given in the first row of Tab. [TJ 

For the other two interactions, in analogy to the YNG interaction [29l [30] , we parametrize each radial 
term of the force ([3T|) as a sum of three gaussians, specifically. 



VANir) = E «^r"exp[-(r//3/'^)2] , (33) 

i=l,3 

for the scalar and spin terms, and 

VUr)^ E wyeM-ir/Plf] , (34) 

1=1,3 

for the tensor term. 

The interaction LNS and the finite range interaction LND are used in calculations where only the 
direct terms of the interaction matrix elements, see Eq. (jl6p . are considered, while we consider both 
direct and exchange terms when we use the finite range interaction we call LNDE. The values of the 
parameters of each interaction have been chosen to reproduce at best the energies of the first 0~, 2~ and 
1~ excited states of . These parameters values are given in Tab. [T]and the momentum dependence 
of the three interactions is shown in Fig. [TJ In Tab. [5J we compare the results of our calculations in 
j^O with the experimental values presented in Refs. [2I1[311[33'. In the table all the energies have been 
renormalized to the energy of the lowest state, which, in our calculations, is always the 0~. 

We have also calculated the spectrum of JfO with the YNG interaction of Refs. [29l [30]. This force 
reproduces the G-matrix interaction built on a Nijmegen A -nucleon potential (34j . The spectrum we 
have obtained with this interaction, is rather different from the experimental one. The lowest energy 
state is a 1+ followed by a 1~ at 0.49 MeV, a 2~ at 13.34 MeV. The first 0" state appears at 17.94 
MeV. We expected such a bad result, since our theory requires interactions which effectively consider 
also many-body effects not included in the G-matrix calculations [35j . For this reason a straightforward 
comparison between our effective interactions with more microscopic interactions is not appropriate. 



4 Results 

Before presenting the results obtained by applying our model to a^C , JfO , j^Ca , A'Zr and ^"^Pb hyper- 
nuclei, we make a comment about the presence of some ground state correlations in our model. In our 
approach, these correlations are taken into account by the Z amplitudes. These effects are related to the 
RPA- like description of the ground state [15], and they describe a limited set of correlations. For example 
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Figure 2: Differences between RPA-like and TDA energy eigenvalues of a'O and J^^Vh hypernuclei calculated 
with the LN(5 interaction. 



they do not describe the short-range correlations [36], and they do not consider the partial occupation 
probability of single particle states both below and above the Fermi surface. We have already pointed 
out that in our model nucleonic correlations are not considered. 

A quantitative indication of the relevance of these RPA-like ground state correlations, can be obtained 
by comparing our results with those of the TDA theory obtained by solving Eq. ([7]) without the term 
which multiplies the Z amplitude. Since we work with a discrete configuration space, the number of 
eigenstates corresponding to positive energies obtained in RPA and TDA is the same. For this reason we 
can immediately identify the energy difference of each eigenstate. As an example of the results we have 
obtained, we show in Fig. [2] these energy differences for the and ^°*Pb hypernuclei as a function 
of the RPA energy. These results have been obtained by using the LN^ interaction. In the left panel 
of the figure we show 273 energy differences corresponding to all the solutions we found for all the 
positive and negative parity multipole excitations up to J = 4. In the case of ^°*Pb we show 5024 energies 
differences corresponding to multipole excitations up to J = 10. 

We observe that only in few cases the energy differences are larger than 10 keV. Specifically, for 
the states of our interest, below 10 MeV for and below 4 MeV for ^"^Pb , these differences are even 
smaller. The results of our RPA-like theory and those of the TDA are quantitatively equivalent. However, 
since the RPA is the complete and consistent theory describing linear combination of one-particle, one- 
hole excitations |15j . we used our RPA-like model to obtain all the results presented in the remaining 
part of the article. 

We start the presentation of our results by showing in Tab. [2] the spectrum oi . In this table 
also the MF results and the experimental values of Refs. [32 [311 [33] are shown. The energy differences 
between the 0", 2^ and \2 states have been used to adjust the parameters of the A-nucleon forces. 
These parameterizations are not unique. We have found sets of parameters values which produce very 
collective 0+ states at low energy, in some case this was even the lowest energy state. While there are 
no experimental evidences of such low energy 0+ state, there are indications of a 0+ state around 10-11 
MeV [371 [33] . For this reason we have used forces predicting a 0+ state around 1 1 MeV. The energies of 
the other states are a prediction of our model. 

The MF energies are already rather close to the experimental values and our RPA-like results change 
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Figure 3: Transition form factors, Eq. (|35p . obtained in MF approximation and calculated with the three residual 
interactions used in this work, as a function of the momentum q of the operators (|21|) and (|24p . For the unnatural 
parity excitations 0~ and 2~, the operator inducing the transition is indicated, see Eq. (I24p . The subindexes 1 
and 2 of the 1~ states indicate the lower and higher energies states respectively, see Tab. [S] 



them by about 4% on average. The experimental spHtting between the two 1~ states is not well repro- 
duced. Better agreement could be obtained with a rearrangement of the force parameters, or by including 
additional force terms such as the spin-orbit or velocity dependent terms. 

The structure of the 0+ state is somehow more collective than those of the other states of the spectrum. 
This feature becomes more evident by observing the transition matrix elements from the ground state to 
the hypernuclear excited states. In Fig. [3] we show the transition form factors 

F{q) = |< ^o\T{q)\^u >|' , (35) 

for the negative parity states of the spectrum and for the 0+ and the 2+ states. The transition operators 
are defined in Eqs. (|2ip and (|24p for natural and unnatural parity states respectively. For unnatural 
parity excitations, we use transition operators characterized by L = J ± 1. In the case of the 0~ state 
only the Tj+i is active, while for the 2^, both operators are active. In this last case, however, we found 
that the Tj+i form factor is zero in mean-field approximation, and four orders of magnitudes smaller 
than the Tj_i form factor when the residual interaction is active. For this reason we neglect it in our 
discussion. In Fig. [3] the transition form factors for the two 1~ states are labeled with the subindexes 1 
and 2 to indicate the lower and higher energy states respectively, see Tab. [2] Finally, in the lower panels 
we show the transition form factors for the 2+ and 0^ states. 

The behavior of these form factors as a function of q is rather interesting. Only the 0"*" transition 
shows a maximum at values of q near zero. For this reason substitutional experiments such as {K^ , tt") 
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Table 3: Electromagnetic transitions between various states of a'O , in Weisskopf units, compared to the shell 
model results of Ref . [39] . 

are better suited to investigate the 0^ states than experiments virhich produce hypernuclci with a rather 
large value of the momentum transfer, such as (e,e'if+) [50] . 

In Fig. [3] the MF transition form factors are compared with those obtained in our RPA-like calcula- 
tions with the three interactions. For the cases we have investigated, we have found that the presence of 
the interaction does not modify sensitively the MF results, but for the O"*" state. In this last case, while the 
LNi5 result is similar to the MF one, the two finite-range interactions produce rather different form fac- 
tors. In the MF model the particle-hole excitation generating the lowest O"*" state is the (1pi/2)a(1pi/2)^^ 
transition. This is also the main transition (W^=0.998) obtained in the calculation with the LNS inter- 
action. The situation changes slightly when the LND interaction is used, since a relevant contribution of 
the (l?33/2)A(lP3/2)n"^ transition (W — 0.121) is associated to the still dominant (1pi/2)a{1Pi/2)^^ tran- 
sition (ly =0.987). The calculation with the LNDE interaction produces a 0+ state with a rather different 
structure. The main transition is now the (lp3/2)A(l?33/2)n ^ (14^=0.843) with a large contribution of the 
(lPi/2)A(lpi/2)ra"^ component (1^=0.366) and also of the (1si/2)a(1si/2)^"^ component (ly =-0.327). 

We obtain another indication of the relative collectivity of the O"*" state by studying the electromagnetic 
transitions between the excited stated of . The B-values of these transitions are shown in Tab. [3l 
The transitions forbidden in MF calculations have relatively small values when the residual interaction 
is switched on. All the Ml transitions are modified on the second significant figure by the residual 
interaction. We use three digits to show the difference between the results obtained with the various 
interactions. The situation changes when the 0"*" state is present. The B-values of the (O"*" lj~) 
transition are quite different in the various calculations. The LNS result is similar to that obtained in 
MF approximation, while the values obtained with the LND and LNDE interactions are 38% and 48% 
larger than the MF value. This result is again a consequence of the collective structure of the 0"*" state. 

In Tab. Owe compare some of our results with those obtained in shell model calculations [35] • The 
order of magnitude of the B-values is the same in the two calculations, but the differences are remarkable. 
Our results are roughly a factor two smaller for the 1J~ — *■ 0~ transition and a factor two larger for the 
I2 li transition. The difference between the results for the other transition is 17%. It would be 
interesting to compare with experimental values. Unfortunately, while the experimental identification 
of the energies of the spectrum is quite accurate ^40}, the measure of the B-values of the transition 
requires higher statistics. The technological improvements are very promising and these measurements 
will become feasible in the near future. 

In Tab. |4]we compare the spectrum with the energies measured in the E369 [4] and FINUDA 
[6] experiments. The uncertainty in the energy resolution of the E369 experiment is about 1.5 MeV, 
and that of the FINUDA experiment is 1.29 MeV. The angular momentum and parity assignment to the 
energies measured in the FINUDA experiment is our guess. 

The presence in the experimental spectrum of two 1^ states, that we have called 1^3, is not predicted 
by our model. Shell model calculations [H] [42] indicate that these states are produced by neutron- 
A transitions where the neutron is in the lpi/2 state, above the Fermi surface in the MF model. The 
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ph 


MF 


LN(5 


LND 


LNDE 


E369 


FINUDA 




(1si/2)a(1?33/2),T^ 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 




(1Si/2)a(1P3/2),T^ 


0.0 


0.26 


0.24 


0.27 


















2.63 


2.54 














6.09 


5.04 


0+ 


(ip3/2)a(ip3/2);:^ 


10.03 


9.61 


8.39 


8.28 


8.12 


7.14 


2+ 


(1P3/2)a(1P3/2);:^ 


10.03 


9.38 


9.33 


9.79 




9.34 


1+ 


{ip3/2)a(ip3/2);:^ 


10.03 


9.69 


9.58 


9.94 


11.00 


11.21 


3+ 


{1P3/2)a(1P3/2);:^ 


10.03 


9.87 


9.80 


9.95 







Table 4: Energies in MeV of the ^^C spectrum calculated in MF approximation and with the three interactions 
used in our work. The experimental values of E369 [3] are those reported in [TT]. The values of the FINUDA 
experiment [6] have been rescaled with respect to the lowest energy value, and the angular momentum and parity 
assignment is our guess. The uncertainty in the energy resolution of the E369 experiment is 1.45 MeV, and that 
of the FINUDA experiment of 1.29 MeV. The MF particle-hole transitions are also shown. 

ground state correlations of our approach do not consider partial occupation probability of the single 
particle levels, therefore we are unable to predict these excitations. A description of ^^C as closed shell 
nucleus is reasonable in kinematic regions where the detailed shell structure is not relevant, such as giant 
resonances or the quasi-elastic region, but for the case we are treating it is inadequate. 
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MF 


LN,5 


LND 


LNDE 


43J 


Ml (2- 


-ID 


0.0 


0.082 


0.082 


0.082 




El (0+ - 


-ID 


0.278 


0.246 


0.421 


0.509 


0.475 


El (2+ - 


-ID 


0.139 


0.135 


0.121 


0.158 


0.133 


El (2+ - 


-2-) 


0.139 


0.170 


0.194 


0.133 


0.094 



Table 5: Electromagnetic transitions between various states of ^^C , in Weisskopf units, compared with the shell 
model results of Ref . [33] . 

Our model predicts the presence of excited states between 8 and 1 1 MeV. These states are all coming 
from the (lp3/2)A(lp3/2)n ^ excitation, and in MF approximation are degenerate at 10.03 MeV. The 
residual interaction removes this degeneracy, and, in general, the excitation energy of each state is lowered 
with respect to the MF value. The 0+ and 2+ states are the lowest ones. The LND and LNDE interactions 
give the 0+ around 8.5 MeV, where E369 finds a peak. The identification of this experimental peak with 
our 0^ is however quite dubious, since the experiment produces hypernuclei by using a probe that transfers 
high momentum values to the system, while the 0^ response has its maximum at small q values. This is 
shown in the lowest right panel of Fig. 01 where the transition form factor of this 0+ state is shown. 

In our search for collective effects we have studied the transition form factors of the various excitations, 
and, in Fig. [4l we show some of them. This figure is analogous to Fig. [Sj but for ^^G . We have selected 
the transition form factors of the positive states degenerate in MF, and the lowest 1" state. For the 3+ 
state we show only the Tj-i form factor, since the Tj+i is zero in the MF model. 

The figure shows that the 2+ and 3+ transitions are essentially unaffected by the effective interac- 
tion, as it is the ground state 1~ transition. The 0+ state, and, much more, the 1+ state, are instead 
rather sensitive to the presence of the residual interaction. The structure of the 1+ state is not collec- 
tive. In our calculation this state is described by the combination of two particle-hole transitions, the 
(1P3/2)a(1P3/2)^^ transitions, which is the dominant one, and the (1pi/2)a(1P3/2)^^ transition. While 
in MF approximation this state is described only by the first transition, the relative weight of the two 
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particle-hole components changes with the residual interaction. The values of the W amplitudes vary, 
respectively, from 0.99 and 0.10 for the LNDE interaction to 0.94 and 0.33 for LNS. These different 
couplings between the two states are sufficient to produce the spreading of the form factors shown in Fig. 

SI 

Also the RPA-like structure of the 0+ state is mostly described by the combination of the (1p3/2)a(1P3/2) 
and {lsi/2)A(lsi/2)n ^ transitions. These are the only two relevant transitions in the formation of the 0^ 
state. The relative weights of these two transitions change by changing the residual interaction, and this 
produces the spreading of the form factors shown in Fig. |4l 

The B-values of some of the 7 transitions for the hypernucleus are shown in Tab. [51 and they 
are compared with the shell model results of Ref. [13]. The 2~ ^ 1^ Ml transition is not present in the 
MF model, since the two states are degenerate. Our RPA-like calculations remove this degeneracy and 
produce B-values which are almost insensitive to the different residual interactions. 

The B-values of the El transitions 2+ 1^ and 2+ ^ 2~ obtained in the MF model are identical, 
since the wave functions of the initial and final states are the same, as it can be deduced from the MF 
transitions shown in Tab. ID The B-values obtained in our RPA-like calculations are rather similar to the 
MF values. The situation changes for the 0+ ^ Ij^ El transition, where our RPA-like results modify the 
MF value by 10, 70 and 100% when the LNS, LND and LNDE interactions are used, respectively. As for 
]yO , also in this case the structure of the 0+ state is extremely sensitive to the residual interaction. In 
our model this state is described by a combination of the (lp3/2)A(lp3/2)n ^ transition, the dominant one, 
and the (lsi/2)A(lsi/2)n "'^ transition. The coefficients of these two particle-hole transitions are changed 
by the residual interaction and this produces the changes in the B-values. The agreement between our 
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results with those of the shell model [33] is more satisfactory than in the case of the j^O hypernucleus. 

In Tab. [6] we show the results relative to the spectrum of the '^Ca. hypernucleus, for energies smaller 
than 10 MeV. The changes of the RPA-like energies, with respect to those of the MF calculations, are 
3.5%, on average. The only exception is the 1^ state which shows changes around 20 %. We have analyzed 
the transition form factors of this state, but we did not find any remarkable effect which distinguishes it 
from the MF result. 



J"" 


ph 


MF 


LN(5 


LND 


LNDE 


m 




(lSl/2)A(lrf3/2),;^ 


0.0 


0.0 


0.0 


0.0 


0.22 




(lSl/2)A(ld3/2),^^ 


0.0 


0.16 


0.21 


0.19 


0.0 


0+ 


(1si/2)a(2si/2)-i 


2.60 


2.50 


2.74 


2.42 


2.92 


It 


(1Si/2)a(2Si/2),T^ 


2.60 


3.07 


3.13 


2.93 


3.44 


2+ 


(lSl/2)A(lrf5/2);;^ 


6.00 


6.04 


6.08 


6.09 


6.19 


3+ 


(lSl/2)A(lrf5/2);;^ 


6.00 


6.10 


6.10 


6.11 


6.01 


2r 


(lP3/2)A(lrf3/2);:^ 


7.81 


7.72 


7.69 


7.77 


7.98 


0^ 


(lP3/2)A(lrf3/2);:^ 


7.81 


7.77 


8.02 


7.07 


9.54 


ir 


(lP3/2)A(ld3/2);:^ 


7.81 


7.95 


7.91 


8.05 


8.78 


3- 


(lP3/2)A(ld3/2);:^ 


7.81 


8.04 


8.24 


7.99 


8.55 


22^ 


(lPl/2)A(ld3/2);:^ 


8.58 


8.61 


8.59 


8.59 


9.53 


^2 


(ipi/2)A(id3/2);:^ 


8.58 


8.74 


8.51 


8.36 


9.27 



Table 6: Energies in MeV of the A'Ca spectrum calculated in MF approximation and with the three interactions 
used in our work. Our results are compared with the shell model results of Ref. [H]. We have considered all the 
states with energy smaller than 10 MeV. The MF particle-hole transitions are shown. 

We observe in Tab. [5] that, in our calculations, the ground state of the j^Ca system is produced by 
the (1si/2)a(1'^3/2)^^ transition which in MF approximation gives two degenerate states characterized 
by angular momenta and parity 1+ and 2^. All our RPA-like calculations remove the degeneracy by 
lowering the 1+ state, while it is the opposite in the calculation of Ref. [44]. In general, the comparison 
between our results and those of Ref. |44l shows an agreement in terms of clustering of states around the 
MF single-particle transitions. However, within the set of states produced by the same MF transition, 
there are some inversions of the multipoles. This happens for the MF degenerate 2^ and 3+ states, 
as well as for the set of degenerate states produced by the (lp3/2)A(ld3/2)n ^ transition and also for the 
(lpi/2)A(lrf3/2)n ^ transition. A measure of the experimental spectrum of j^Ca. would clarify the situation. 



-Zr 


P 


MF 


LNS 


LND 


LNDE 


exp [4] 


(1Si/2)a 


0.0 


0.003 ± 0.003 


0.005 ± 0.009 


0.007 ± 0.012 


0.0 


(1J5i/2)a 


4.994 


5.012 ± 0.061 


5.010 ± 0.058 


4.999 ± 0.031 


4.9 


(lt^5/2)A 


10.749 


10.783 ± 0.077 


10.766 ± 0.060 


10.743 ± 0.038 


11.5 


(1/7/2)a 


16.994 


17.016 ± 0.035 


17.013 ± 0.039 


17.002 ± 0.029 


19.0 



Table 7: Energies, in MeV, of the various states coupled to the lgg/2 neutron hole in A'Zr , calculated in 
MF and in our RPA-like approach for the three A-nucleon interactions used in this work. The energies of the 
RPA-like results are the averages of the various excitation energies of each multipole excitation compatible with 
the indicated single particle transition. The experimental energies are those of the hypernucleus, and the 
uncertainty in the determination of the energy is roughly 1 MeV. 

We have studied the spectra of the ^"Zr and j^^Ph hypernuclei and we have observed that also in 
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these cases the changes with respect to the MF results are rather small, about 100 keV, of the order of 
the 2, maximum 3% on average. 

As an example of our results, we show in Tab. [7]the energies of the excited states of ^°Zr obtained by 
considering the transitions of the I59/2 neutron hole to the A single particle states. This MF picture has 
been used in the analysis of the (tt"*", fsr+) reaction on ^^Y target In the same spirit, we have studied 
the excitation of this neutron hole into the s,p,d, f states of the A with = l\ + l/2. We have considered 
all the possible multipole excitations compatible with a specific neutron-A transition. Obviously, in the 
MF calculation all these states have the same excitation energy. Our RPA-like calculations, however, 
give different energy values for each multipole excitation. In Tab. [7] we compare the MF energies 
for each transition with our RPA-like energies obtained by averaging the results of the various multipole 
excitations, and we also show their standard deviation. We observe that, within the statistical uncertainty, 
RPA-likc and MF results coincide. Also the agreement with the experimental values is good, especially 
if we consider that the data are taken on a ^^Y target, and they have an uncertainty of 1.65±0.10 MeV. 



p 


MF 


LN(5 


LND 


LNDE 


expl 


exp2 


(1si/2)a 


1.633 


1.630 ± 0.001 


1.629 ± 0.002 


1.623 ± 0.008 


1.6 


-3.0 


(1P1/2)a 


4.621 


4.629 ± 0.024 


4.626 ± 0.025 


4.616 ± 0.016 


6.8 


4.6 


(lt^5/2)A 


8.220 


8.237 ± 0.036 


8.230 ± 0.030 


8.216 ± 0.017 


12.84 


8.24 


(1/7/2)a 


12.344 


12.354± 0.020 


12.345 ± 0.010 


12.328 ± 0.014 


18.29 


13.69 


(139/2)a 


16.926 


16.948 ± 0.034 


16.959 ± 0.040 


16.926 ± 0.108 


22.10 


17.49 



Table 8: Energies, in MeV, of the various states coupled to the lii3/2 neutron hole in ^°*Pb , calculated in 
MF and in our RPA-like approach for the three A-nucleon interactions used in our work. The energies of the 
RPA-like results are the averages of the various excitation energies of each multipole excitation compatible with 
the indicated single particle transition. The values shown in the columns expl and exp2 are the experimental 
energies of Ref. [2], rescaled in two different manners (see text). 

The results of an analogous study for the 1*13/2 neutron hole in ^"^Pb coupled with the s,p,d,f,g 
states of the A with ja = ^a + 1/2, are shown in Tab. [H] The RPA-like energies are all compatible with the 
MF ones, within the statistical uncertainty. In the same table we compare our results with the energies 
given in Ref. [2] where the (tt"*", if +) scattering data on a ^°^Pb target have been analyzed in terms of 
transitions of the lzi3/2 neutron hole. Since, in our calculations, the lowest energy of ^"^Pb is produced 
by the transition (1si/2)a(3pi/2)^^, we have positioned the first experimental energy of Ref. [2] at 1.63 
MeV, which is the difference between the 3pi/2 and l?i3/2 neutron single particle energies. The other 
energies have been consistently rescaled, and their values are shown in Tab. [8] by the column labeled 
expl. The agreement is not satisfactory even if we consider that the experimental energy resolution is of 
2.2 ± 0.1 MeV. We then set the position of the second experimental energy at 4.6 MeV, to match our 
second state, and we rescale the other energies. The corresponding values are shown in the column exp2. 
In this case the agreement with our results noticeably improves. The source of disagreement with our 
results is in the energy difference between the first two states. 

Up to now, we have presented our results by considering individually each hypernucleus. To give a 
more general view, we show in Fig. [5] the density of states of the various hypernuclei as a function of the 
excitation energy. The different lines indicate the results obtained in MF and RPA-like approaches by 
using the various residual interactions. In our model, the number of excited states for a given multipolarity 
depends only on the configuration space, and not on the residual interaction. For this reason, differences 
between the lines of each panel of the figure indicate the relevance of the A -nucleon residual interaction 
which shifts the excitation energies with respect to the MF results. For each hypernucleus, we have taken 
into account all the possible multipole excitations which provide a solution below 20 MeV. 

A first remark on the results shown in Fig. [5] is that the density of states in ^"^Pb is almost ten 
times larger than that of ^■^C or ^^O . This is an indication of the energy resolution required to do 7 
spectroscopy in heavy hypernuclei. A second remark, more related to the aim of our investigation, is that 
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the differences between tfie various lines seem to be more evident for light hypernuclei than for the heavier 
ones. This could indicate that the MF approach becomes more reliable the heavier is the hypernucleus, 
or in other words, that collective effects are more inhibited in heavy hypernuclei than in the light ones. 

In order to investigate this issue, we have calculated the normalized energy shift with respect to the 
MF energy. 



AE 




J,RPA 



J,MF 



J,RPA 



J,MF 



(36) 



where we have indicated with Nj the number of multipole excitations we have considered, and the index 
i runs on all the S solutions having the same multipolarity. The above equation is meaningful since, 
as we have already mentioned, our model produces the same number of excited states in both MF and 
RPA-like calculations, therefore we can identify the energy shift of each MF solution. We show in Tab. [H 
the AE values obtained for the various A-nucleon interactions. These results indicate that the deviations 
from the MF model are less important, the heavier is the nucleus. 

We have compared the values of the hyperonic energy differences of Tab. [9] with analogous values 
calculated for excitations of ordinary nuclei. The RPA calculations for the doubly magic nuclei have 
been done with the zero-range Landau-Migdal interaction of Rcf . 45j . This comparison indicates that 
collectivity in hypernuclei is much smaller than in ordinary nuclei. The AE values for ordinary nuclei 
are about five times larger than the analogous values obtained for hypernuclei. 



16 



J. 



max 



LN6 



LND 



LNDE 




5 
4 
4 
6 

10 



0.23 ± 0.07 
0.33 ± 0.10 
0.15 ± 0.04 
0.10 ± 0.02 
0.10 ± 0.02 



0.26 ± 0.08 
0.30 ± 0.09 
0.12 ± 0.04 
0.13 ± 0.03 
0.13 ± 0.02 



0.22 ± 0.06 
0.11 ± 0.03 
0.24 ± 0.07 
0.09 ± 0.02 
0.09 ± 0.01 



Table 9: Renormalized energy differences Eq. (|36|l . We have also indicated the maximum value of the angular 
momentum of the multipole excitation. 

5 Summary and conclusions 

In this article we have investigated the possible presence of collective phenomena, beyond the MF model 
description of hypernuclei. To this purpose we have developed a particle-hole model where the formation 
of the hypernucleus is described as the excitation of a nucleon into a single particle state describing 
the A , and the hypernucleus state is a linear combination of elementary single particle excitations, the 
coefficients of the combinations being determined by the residual A -nucleon interaction. Our model is 
inspired by the RPA theory and considers some ground state correlations which can be described in terms 
of A -nucleon excitations. Purely nucleonic correlations, such as those considered by the traditional RPA 
theory, are not included. 

Our calculations have been done by using a discrete single particle configuration space, and with three 
A -nucleon interactions with different properties. Since we have been interested in the low energy part of 
the spectrum, we have checked the stability of our results with respect to the size of the configuration 
space, and therefore the scarce relevance of a detailed treatment of the continuum. 

When possible we have compared our results with those obtained by other calculations. We obtained a 
reasonable agreement in and ]f O , not only for the energies [HI 112], but also for the electromagnetic 
B-values with shell model results [331 [33] • For the j^Ga. hypernucleus, the comparison with the shell 
model spectrum of Ref. 04] is not so satisfactory. We have the same clustering of states in the spectrum, 
but we often obtain an inversion of states within the cluster. 

Concerning the comparison with the experimental data, we obtain a reasonable agreement with the 
energies measured in ^^O . Our model is unable to predict two states observed in because it does 
not consider the partial occupation of the lpi/2 neutron state. The energy resolution of ^°Ca data [371 [T] 
does not allow a meaningful comparison with our results. For the ^°Zr and j^^Ph hypernuclei we made 
a comparison with data analyzed in terms of excitation of a single neutron [4] [2] . The agreement with 
^°Zr data is satisfactory, while for the j^^Ph case there are some problems regarding the spacing between 
the first two states. 

To investigate the presence of collective effects in the spectra of hypernuclei, we have calculated the 
form factors of the transition producing the hypernucleus and we have compared MF and our RPA-like 
results. In few cases we have found remarkable differences, however, in all the cases but one, these 
differences could be explained by considering the excited state as combination of only two single particle 
transitions. This is not really what we would call a collective behavior of the system. The only situation 
standing out of this trend is that of the 0+ excitation in \^0 around 10 MeV, which seems to have a 
more collective structure. 

The energy shifts with respect to the MF results are indicators of effects coming from the A -nucleon 
interaction, and possibly of collective phenomena. Densities of states and average energy shifts do not 
show remarkable features beyond the MF approach. In our calculations the differences between MF 
and our RPA-like results become smaller the heavier is the hypernucleus. This is an indication that 
perturbation effects produced by the A on the nuclear systems are better averaged in a heavy system. 

The average energy shifts for hypernuclear excitations are, at most, 20% of the shifts obtained for 
the excitation of a purely nucleonic system. This result is almost independent of the number of nucleons 
involved, and indicates that in hypernuclear excitation, effects beyond the MF are much less important 
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than in the excitation of the ordinary nuclear system. 

Our model can be further improved in several ways. One of them consists in performing self-consistent 
calculations, where the single particle basis is generated in a Hartree-Fock approach, with the same A- 
nucleon interaction used in the RPA-like calculations. Other improvements can be obtained by considering 
nucleonic correlations in the ground state. This could be done in terms of hole-particle excitations, as in 
the ordinary RPA, in terms of partial occupation probability of the nucleonic single particle states, as it 
happens in the quasi-particle RPA, and also by considering explicitly the short-range correlations as in 
Ref. IMI- 
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Appendix 

In this appendix we give the expressions of the matrix elements and ([TT]) of the secular RPA-like 
equations for all the terms of the interaction p^ . 
For the scalar term V^^y obtain the expressions 



J, dir 
^c,AhA'h' 



dr2rlRl,{r2)RA'{r2)jj{qr2) 



(_lpA+iA'+l , 



47r 



-]hJh]A'Jh' 



3A J 
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2 



Jh 

-1 



£.{lA + lh + J)£.{Ia' +lh' + J) 



Jh' 
1 

2 



J JA' 

-i 



(37) 



and 



J, exch 
'^cAhA'h' 



dqq^^Aiq) J dnrlRl{n)RA' {nMqn) 



dr2 r2Rh'{r2) Rh (?'2 )ji {qr2 ) 



-jAjhjA'jh'l 



3A jh 
jh' JA' 



JA I JA' 



J 



I > ^{Ia + ia' + m{ih + ih' + 1) 



Jh' 

1 

2 



Jh 

1 

2 



(38) 



In the above equations we have indicated with R{r) the radial part of the single particle wave functions, 
with ji{qr) the spherical Bessel functions, and with the traditional symbols the Wigner three and six j 
coefficients [T51. 
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The analogous expressions for the spin-dependent terms are 
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For the tensor components of the interaction we obtain: 
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for the direct term, and for the exchange term we get the expression 
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dr2 rlRl,{r2) Rh {r2 )3h (9^2 ) 

lAihiA'ih' jAjhjA' jh'lili 
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h h 2 \ / /a ^1 

Q Q )\ Q 

h h 2 

1 1 L 








JK Jh 

jh' jh' 



Ih' 


1 

! 


jh 


Ia' 


2 


jA 


k 


1 


L 



J 

L 



(42) 
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